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Abstract. Inverse scattering methods capable of compressive imaging are proposed and ana- 
lyzed. The methods employ randomly and repeatedly (multiple-shot) the single-input-single-output 
(SISO) measurements in which the probe frequencies, the incident and the sampling directions are 
related in a precise way and are capable of recovering exactly scatterers of sufficiently low sparsity. 

For point targets, various sampling techniques are proposed to transform the scattering matrix 
into the random Fourier matrix. Two schemes are particularly interesting: The first one employs 
multiple frequencies with the sampling angle always in the back-scattering direction resembling the 
synthetic aperture (SA) imaging; the second employs only single frequency with the sampling angle 
in the (nearly) forward scattering direction in the high frequency limit, resembling the setting of 
X-ray tomography. 

The results for point targets are then extended to the case of localized extended targets by 
interpolating from grid points. In particular, an explicit error bound is derived for the piece-wise 
constant interpolation which is shown to be a practical way of discretizing localized extended targets 
and enabling the compressed sensing techniques. 

For distributed extended targets, the Littlewood-Paley basis is used in analysis. A specially 
designed sampling scheme then transforms the scattering matrix into a block-diagonal matrix with 
each block being the random Fourier matrix corresponding to one of the multiple dyadic scales of the 
extended target. In other words by the Littlewood-Paley basis and the proposed sampling scheme 
the different dyadic scales of the target are decoupled and therefore can be reconstructed scale-by- 
scale by the proposed method. Moreover, with probes of any single frequency to the coefficients in 
the Littlewood-Paley expansion for scales up to uj/(2-k) can be exactly recovered. 



1. Introduction 
Consider the scattering of the incident plane wave 

(1) ^(r) = e iwr - d 

by the variable refractive index n 2 (r) = 1 + u(r) where d is the incident direction. The scattered 
field satisfies the Lippmann-Schwinger equation [9j 

(2) u s {r) = oj 2 j v(r') (V (r') + u s (r')) G(r, r', u)dr', reR d ,d = 2,3 

where G(r, r' ,ui) is the Green function of the operator —(A + u; 2 ). We assume that the wave speed 
is unity and hence the frequency equals the wavenumber uj. 
The scattered field has the far-field asymptotic 

iu>\r\ 

(3) n s (r)= /2 (A(r,d,o;) + O(|r|- 1 )), f = r/|r|, 
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Figure 1. Far- field imaging of discrete or extended scatterers 



where A is the scattering amplitude. In inverse scattering theory, the scattering amplitude is the 
measurement data determined by the formula [9] 



The main objective of inverse scattering then is to reconstruct the medium inhomogeneities v 
from the knowledge of the scattering amplitude. In Part I [12J and this paper the target to be 
imaged consists of a finite number of "point scatterers. And the main techniques for reconstruction 
are from theory of compressed sensing. In |12j we analyze the one, but high, frequency imaging 
method with the single-input-multiple-output (SIMO) and multiple-input-multiple-output (MIMO) 
measurements in which for every incident plane wave the scattering amplitude is sampled at multiple 
directions independent of the incident wave. 

In this paper the focus is on the multi-shot single-input-single-output (SISO) measurement in 
which for every randomly selected incident plane wave the scattering amplitude is sampled at only 
one direction correlated with the incident wave. 

Our motivation for this alternative imaging method is practical as well as theoretical. On the 
theoretical aspect, the analysis of the high-frequency SIMO/MIMO method employs the coher- 
ence theory of compressed sensing which deals with only random targets under a suitable sparsity 
constraint. On the other hand, the multi-shot SISO method proposed in this paper is amenable 
to the restricted isometry theory of the random Fourier matrix which guarantees reconstruction 
for all targets under the weakest known sparsity constraint. On the practical aspect, the present 
method can achieve a comparable performance with a much lower frequency or bandwidth (Figure 
[7]). Moreover, the case of extended targets can be treated by either interpolating from grid points 
or using the wavelet basis in this approach. The main drawback, though, of the present approach, 
in comparison to that of |12j . is that the multiple scattering effect is not accounted for. 

In Sections [2j we discuss the case of point scatterers and propose several sampling schemes to 
transform the scattering matrix into the random Fourier matrix which is amenable to the com- 
pressed sensing techniques. One scheme employs multiple frequencies with the sampling angle 
always in the back-scattering direction resembling the synthetic aperture (SA) imaging; another 
scheme employs only single frequency with the sampling angle in the (nearly) forward scattering 
direction in the high frequency limit, resembling the setting of X-ray tomography. We then extend 
these results to the case of localized extended targets by interpolating from grid points in Section [3} 
In Section [4] we analyze the case of distributed extended targets using the Littlewood-Paley basis 



(4) 
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and propose a sampling scheme to block-diagolize the scattering matrix. Each block is in the form 
of random Fourier matrix and corresponds to one dyadic scale of the target. Hence our method has 
the capability of imaging the target scale-by-scale by the compressed sensing techniques. Moreover, 
the coefficients in the Littlewood-Paley expansion for scales up to lu/(2tt) can be exactly recovered 
by using probes of any single frequency to. We numerically test these sampling methods and com- 
pare their success probabilities in Section [5j We conclude and comment on the issue of resolution 
in Section [HI 



2. Point scatterers 

For the simplicity of notation, we will focus on two dimensions below. 
We consider the medium with point scatterers located in a square lattice 

C = {rj = (xi, zi) :i = 1, m} 

of spacing I. The total number m of grid points in C is a perfect square. Without loss of generality, 
assume Xj = Zj = j 2 i where j = (ji - l)y/m + j 2 and j l7 j 2 = 1, \fm. Let Uj, j = 1, m 
be the strength of the scatterers. Let S = jr^. = (xj. ,z%- ) : j = 1, ...,s} be the locations of the 
scatterers. Hence Vj = 0, Vr^ g" S. 

For the discrete medium the scattering amplitude becomes a finite sum 

2 m 

(5) A(f,d, W ) = ^E^fe) e ^ rj " f - 

j'=i 

Unlike [12] which covers both linear and nonlinear scattering, here we work exclusively under the 
Born approximation in which the exciting field u{tj) is replaced by the incident field v}(rj); unlike 
|12j which deals exclusively with one frequency, here we will work with multiple frequencies. 

Let di,ri,l = 1, n be various incident and sampling directions for the frequencies ui,l = 1, n 
to be determined later. Define the measurement vector Y = {Y{) G C n with 

(6) Yi = ^A(ri,di,u;i), 1 = 1, ...,n. 

The measurement vector is related to the target vector X = (vj) G C m by the sensing matrix $ as 
Y = <&X. Let 0i, 9i be the polar angles of d^r^, respectively. The (/,j)-entry of $ G c nxm is 

(7) g-iwjfjTj-giwjdj-rj _ e tWt£(j 2 (sm0i-Bm§i)+ji(cosOi-cos9i)) ■ _ r-^ _ -|\ _|_ ^ 

Note the all entries of $ have unit modulus. As in [IL2] we reconstruct X as the solution of the 
^-minimization, called Basis Pursuit (BP): 

(8) min||Z||i s.t. <f>Z = Y 

which can be solved by linear program or by various greedy algorithms [2 [81 [T7] . 
In the presence of noise/error E of size e as in 

(9) Y = $X + E, ||£|| 2 <£ 
Q is replaced by the relaxation scheme 

(10) min ||Z||i, s.t. [|y-*Z|| 2 <e. 

A fundamental notion in compressed sensing under which BP yields the unique exact solution 
is the restrictive isometry property (RIP) due to Candes and Tao [6j. Precisely, let the sparsity s 
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of a vector Z S C m be the number of nonzero components of Z and define the restricted isometry 
constant 8 S to be the smallest positive number such that the inequality 

k(1 - S a )\\Z\\l < \\*Z\\l< K (l + 6 a )\\Z\\l 

holds for all Z G C m of sparsity at most s and some constant k > 0. For the target vector X let 
denote the best s-sparse approximation of X in the sense of L 1 -norm, i.e. 

X^ = argmin \\Z - X\\i_, s.t. \\Z\\ < s 

where \\Z\\q denotes the number of nonzero components, called the sparsity, of Z. Clearly, X^ 
consists of the s largest components of X. 

Now we state the fundamental result of the RIP approach [3] which is an improvement of the 
results of 01 El. 

Proposition 1. [3J Suppose the restricted isometry constant of & satisfies the inequality 
(11) 5 2s <V2-l 



with k = 1. Then the solution X by BP (10) satisfies 

(12) H^-^lh < C lS - 1/2 \\X -X^\\x + C 2 e 

for some constants C\ and C 2 . 

Remark 1. For general ft / 1, one can consider the normalized version of |5p 

x =y = ^=3>X + —=E 



£ 



by which it follows that 

(13) \\X-X\\ 2 < ds-VHx - X('\ + C 2 

We wish to write the (I, j)-entry of the sensing matrix in the form 

(14) e w(3ift+i26) > i = (il _i )v ^ + i2) ^,^ = 1,...,^, l = l,...,n 

where Q are independently and uniformly distributed in [—1, 1] in view of the following theorem. 
Proposition 2. [15] Suppose 

(15) ^— > C<T 2 o- In 2 cr In m In-, a €(0,1) 

Inn a 

for a given sparsity a where C is an absolute constant. Then the restricted isometry constant of 

the matrix with entry (14) satisfies 

8, <S 

with k = n and with probability at least 1 — a. 

See [H d [H] for the case when belong to the discrete subset of [—1, 1] of equal spacing 
2/y/m. 



To construct a sensing matrix of the form (14) we proceed as follows. Write (£/, Q) in the polar 
coordinates pi , as 



(16) te,0)=«(cos^,sin^) ; Pl = \J^i+Qi <V2 
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and set 



w,(cos Qi — cos Oi) 
Ui(sm0i - sin §i) 



COS 01 

sin mi 



where is a parameter to be determined later (30). Equivalently we have 

01 — 01 ■ 01 + &i 
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01 + 01 



VLpi cos ( 



£lpi sin i 



(17) — V^wz sin — - — sin ■ 

(18) V2uji sin ^ - ^ cos ■ o 

This set of equations determines the single-input-(#z, w,)-single-output-#z mode of sampling. 

The following implementation of ( 17 )-( 18 ) is natural. Let the sampling angle 0\ be related to 
the incident angle Q\ via 

(19) 0i + 0i = 2fa + 7T, 
and set the frequency w, to be 

ttpi 



(20) 



V2 



sm ■ 



Then the entries of the sensing matrix $ have the form 



(21) 



31,32 = l,...,Vm. 



1, ...,n, 

By the square-symmetry of the problem, it is clear that the relation ( 19 ) can be generalized to 
(22) 0i + 0i = 20, T) G Z. 

On the other hand, the symmetry of the square lattice should not play a significant role and hence 
we expect the result to be insensitive to any fixed r\ G K, independent of as long as (20) holds. 
Indeed this is confirmed by numerical simulations below (Figure [5]). 
Let us focus on three specific measurement schemes. 



Scheme I. This scheme employs ft— band limited probes, i.e. uji G [—0,0]. This and (20) lead to 
the constraint: 



(23) 



sm 



01-01 



> 



V2' 



The simplest way to satisfy (19) and (23) is to set 

(24) <h = e, = e, + 7r, 

(25) un = -j, 

I = 1, ...,n. In this case the scattering amplitude is always sampled in the back-scattering direc- 
tion. This resembles the synthetic aperture imaging which has been previously analyzed under the 
paraxial approximation in |14j . In contrast, the forward scattering direction with 0\ = Q\ almost 



surely violates the constraint ( 23 ) . 

Scheme II. This scheme employs single frequency probes no less than Q: 
(26) ui = 7O, 



7> 1, I 
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1, ...,n. 




Figure 2. Scattering from extended targets 



To satisfy ([22} and ([20j) we set 

(27) e i = (j)i + r ^ + arcsin Pl 



2 7 V2 

(28) 0, = ^ + IE _ arcsin A 

with i] £ 7*. The difference between the incident angle and the sampling angle is 

(29) 01-01 = 2 arcsin ''' 



which diminishes as 7 — > 00. In other words, in the high frequency limit, the sampling angle 
approaches the incident angle. This resembles the setting of the X-ray tomography. 



Scheme III. This scheme employs probes of unlimited frequency band. Let 9i be n arbitrary 
distinct numbers in [— ir, ir] and let Oi and u>i be determined by (22) and (20), respectively. The 
possibility of having a small divisor in ( 20 ) renders the bandwidth unlimited in principle. 
The following result is an immediate consequence of Proposition [T] and Proposition [2] 

Theorem 1. Let £z,0 be independently and uniformly distributed in [—1, 1] and let (pi,(j>i) be the 
polar coordinates of (£/,C/)> i- e - 

(&,Cz) = «( cos 4>h sin <t>l)- 



Let the probe frequencies oji, the incident angles 0[ and the sampling angles 9i satisfy (19) and 



(20), for example, by Scheme I, II or III. 
Suppose 

(30) 



TT 



and suppose (15) holds with a = 2s and any 5 < \[2 — 1. Then (11) with k 
matrix $ and the bound (13) holds true with probability at least 1 — a. 



n is satisfied for the 
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3. Localized extended targets 



In this section we consider extended targets that are localized in space. We represent such 
targets by interpolating from grid points and thus extend the preceding results for point targets 
to localized extended targets. We will treat the case of distributed extended targets that are not 
localized in space in the next section. 

Suppose that the target function z/(r) is continuous and has a compact support. Consider the 
filtered version u„ of v 

(31) „(r, = /,„(,-,>(,')*' 

where <7^(r) = n~ 2 g(r/r]) for any integrable filter function g such that J g(r)dr = 1. Clearly, v„ 
tends to v as r] — >■ 0. 

Next we discretize (31) by replacing the integral by the Riemann sum of step size £. We obtain 

(32) u n/ (r) =£ 2 ^2g v (r-tc l )v(£c L ), JcZ 2 

pel 

which is as smooth as the interpolation element g^ is. 

Since v has a compact support, I is a finite set. For simplicity let / be the square sublattice 



I = {q = (91,92) ■ qi,Q2 = i, Vm} 



of total cardinality m. Let j = (qi — l)y/m + Define the target vector X = (Xj) G C m with 
Xj = Let uj\ and d; be the probe frequencies and directions, respectively, and let r; be 

the sampling directions for I = 1, n. Now we write the data vector Y in the form Q with the 
sensing matrix elements 

< 33 > •« = ^5^))I 9 ' (r '- fq)e '"' (d '" ,r ' dr ' ^(«-^ + * 

_ e iwj£q-(d{-fj) 

and the error term E due to the filtered discretization. For sufficiently small rj, £ we may assume 
\\E\\2 < en 1 / 2 for a given e > (see Remark [l] and below). 

The crucial observation is that the sensing matrix (33) is identical to ([7]) with (xj,Zj) = wq 
and any isotropic filter function g. Therefore Theorem [1| holds verbatim for the case of localized 
extended targets formulated above. 

How small must rj and £ be in order to ensure that \\E\\ 2 < en 1 / 2 ? This can be answered roughly 
as follows. First, by the inequality H-EH2 < ||i£||oo\/w it suffices to have H-EHoc < £• Now consider 
the transformation 7", defined by 

from the space of continuous functions supported on m£] 2 to C n , cf. Q. By definition 

E = Tv — Tv V: e 

and we have 

1 1 U — U„ I) 1 1 1 

\\E 00 < 



2vrmin / ^(77^ (d/ - f;))| 
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For the sampling schemes I and II (with 7 = 1) of Section [2j we can give an explicit bound for 
the discretization error 



■^ 00 < 



1 



2 il^ ~~ ^HlllS HL°°([-27?r2,2 ? ?f2] 2 )- 
where || • \\L°°([-2ri^l,2r)^l} 2 ) denotes the L°°-norm of functions defined on [—2r]Q,2rjfl]' 2 . 

Theorem 2. Consider the sampling schemes I and II (with 7 = 1 ) of Section^ Assume (30) and 

2tts 



(34) 



< 



\9 llL° c ([-2r?Q,2r?C] 2 ) 



Suppose that (15) holds with a = 2s and any 5 < \[2 — 1. Then the bound (IS) with k = n holds 



true with probability at least 1 — a. 

Remark 2. For example, consider the smooth isotropic filter function 

1 

9(t) = -e~ 2 



with g(k) = <?(k). With the choice rj = i, the condition (34) can be formulated as 

(35) Ik-^lli < e^e. 

For another example, consider the indicator function g on the unit square [—5, ^] 2 - The resulting 
interpolation is the piece-wise constant approximation. Then 

k 

' MM 

g(ki,k 2 ) 

Setting again n = t we obtain the condition 

(36) ||f - VijWi < 



2 sin ^ sin \ 



7T kl 



7T 2 e 



2sin 2 (vr/V2) 



which appears to be much more useful than (35). 



4. Distributed extended targets 

In this section, we consider extended targets represented by square-integrable functions v(x,z). 
To this end we use the Littlewood-Paley basis 



(37) 



^n = l (27r)_1 - ^ < 1^1, ICI < 2vr 

\ 0, otherwise 



or 

(38) V( r ) = (vr 2 xz) _1 (sin {2irx) — sin (vrx)) • (sin (2irz) — sin (ttz)). 
Then the following set of functions 

(39) V P ,q(r) = 2-^+P 2 )/V(2- p r - q), p, q G I? 
with 

2 t p r = (2- pi x,2~ P2 z) 



forms an orthonormal wavelet basis in L 2 (M 2 ) [TT]. In terms of the Littlewood-Paley basis we have 
the expansion 

(40) V(X,Z)= ^P.qVW 2 ^)- 

p.qGZ 2 

With the incident fields 

4(r) = e^ r - d *, k = l,...,n 
we have from Q, ^ and (40) that 

(41) y fc = 2vr £ 2^ + ^)/ 2 z, p , q e^ 2P ^- f '=)-^(a; fc 2P(r fe -d fe )), k = l,...,n. 

Define the sensing matrix elements to be 

(42) d> M = -^—^(^(h - d fc ))e^ 2P ( dfe "^- q 

and let = [3>fc,z], where d k , , ui k are given below. 
Let 

pi-l p 2 -l 

Z = (2mj + l) 2 + (ft + m p )(2m p + 1) + (q 2 + m p + 1), |q|oo < m p , |p|oo<P*, 

ji=-p* h=-v* 
p'i-i P2-1 

fe = J] X) (2nj + l) 2 + ((?i+n P 0(2n p , + l) + (^ + Hp' + 1), 

q 1 00 < rip', |p'|oo < P*, 

h=-p* h=-v* 

for some m p ,ra p ,p* £ N, be the column and row indexes, respectively, of It is important to 
keep in mind how k and I are related to (p', q') and (p, q), respectively, in order to understand the 
scheme described below. 

Let Cfc be independent, uniform random variables on [—1, 1] and define 

* / 1 + 6, & e [0,1] 

" Wfc 2fi 1-1 + 6, & €[-1,0] 

; f 1 + a, a e [0, 1] 

u k 2P2 \-i + Cfc, CfcG[-l,0]' 

Suppose a k , (5 k G [—1, 1] for all p', |p'| < p*. This holds true, for example, when the frequencies Uk 
satisfy the constraint 

(45) w fe 2 p 'i > 2vr, oj k 2 p ' 2 > 2vr. 

Let (pt,4>k) be the polar coordinates of (ctk,(3k)- 

As before, let 8k, &k be the angles of incidence and sampling, respectively, which are chosen 
according to 

(46) - 2 sin — - — sin — - — = a k = Pk cos <pk 

(47) 2 sin — - — cos — - — = (3 k = p k sin <p k 



(43) a k 

(44) fa - 



in analogy to (17)-(18). This means 

(48) - u: k 2»{i k - d fc ) = 0J k (2 Pl a k , 2 P * ft) = ir(2^ (sgnfe) + 6), 2 P2 -^(sgn(C fc ) + &))■ 
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By (48) and the definition of ip it is clear that &k,i are zer ° if P 7^ p'- Consequently the sensing 
matrix is the block-diagonal matrix with each block (indexed by p = p') in the form of random 
Fourier matrix 

1 



(49) 

Let X = {Xi) with 
be the target vector. Let 



k,i 



2n D + l 



X l = 2vr(2n p + l)2^^ 2 u pA 



p* 



m 



E E + 

ji=-p* h=-p* 

We can then write the measurement vector Y = $>X where 3? £ C nxm . The above observation 
means that the target structures of different dyadic scales are decoupled and can be determined 
separately by our approach using compressed sensing techniques. 
To solve (46)- (47) we consider 

(50) 

and 



(51) 



9 k + e k = 2<j) k + vr 

Ok — Ok 



2 sin 



Pk- 



To satisfy the relationships (50) and (51) we set 

(52) e k = 

(53) § k = 



. 7T . p k 

k + - + arcsm — 

.1* ■ Pk 

k + - - arcsm — 



analogous to (27) and (28). 



In the case of extreme anisotropy (needle-like structure), say 2 Pl P2 3> 1, (43)- (44) implies 



and hence 
(54) 



+ e k wo 



by (50). On the other hand if 2^"^ < 1 then 



and hence 
(55) 



±7T 



Relations (54) and (55) are reminiscent of Snell's law of reflection if 6k are interpreted as the 
reflection angles. 

Using the RIP for the random Fourier matrix of each block, we obtain the following theorems 
analogous to Theorem [T} 
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Figure 3. Success probabilities for Scheme I with (58). The label indicates the 



value of n. As the recipe (24) is increasingly violated, the performance degrades 
accordingly. 



n p and a 



2s p and any 5 < v2 — 1 . 



Theorem 3. For each p, |p| < p*, let (15) be satisfied forn 
Suppose 

(56) co k > 2vr-2 max ( pi ' P2 \ k = l,...,n p . 

Let be independently and uniformly distributed in [—1,1]. Let the incident and sampling 



angles be determined by (4-3)-(44), (4G)-(4 r l)- 



Then (11) is satisfied for the p-block of the sensing matrix with probability at least 1 — a. Con- 



sequently the solution X p by BP (10) satisfies 



(57) \\X P -X } 



pl|2 



< 



Cis p 1/2 ||X F 



^ p) ||i + C 2 e 



for the same constants C\,C2 as in (13) 



Remark 3. The condition (56) implies that the wavelengths are no larger than the scales under 
interrogation, consistent with the classical resolution criterion. 



Theorem^ allows reconstruction with probes of single sufficiently high frequency 



UJ 



# 



The beauty of the theorem, however, lies in the fact that with probes of any frequency oj the coef- 
ficients in the Littlewood-Paley expansion for scales up to lo/(2it) can be recovered (exactly in the 
absence of noise). 

5. Numerical results 

Greedy algorithms have significantly lower computational complexity than linear programming 
and have provable performance under various conditions. For example under the condition ^3 S < 
0.06 the Subspace Pursuit (SP) algorithm is guaranteed to exactly recover X via a finite number 
of iterations [10 1 . We have used SP for reconstruction in all our simulations with the following 

li 




Figure 4. Success probabilities for Scheme II (a): rj = l,fj = 1, 1/2, 1/4 with 7 = 1, 20. 




Figure 5. Success probabilities for Schemes II (b): rj = fj = 1, 1/2, 1/4 with 7 = 1, 20. 



parameters: m = 2500, £ = 1,Q = 7r/v2,n = 100. The probability of recovery is calculated by 
using 1000 independent runs. In Figures [3]{7j the vertical axis is for the probability of recovery and 
the horizontal axis is for the number of point scatterers. 



To test Scheme I numerically, we use (25) and 



(58) 0j = 0i + »77r, 77 = 1,1/2,1/4,1/8 

to see if the deviations from (24) have any impact on performance. Their probabilities of recovery 
are plotted as a function of the sparsity in Figure [3j Clearly, the performance deteriorates rapidly 

12 





Figure 7. Solid curves are the success probabilities for the SIMO measurement at 
7 = 1, 20, 200 and the dashed curve is the SISO Scheme II at 7 = 1. 



as the difference between the sampling and incident angles decreases. This is due to increasingly 



more frequent and more severe violation of (23) as a result. In other words, the backward scattering 



direction is the optimal sampling direction for Scheme I. 



Likewise, to test Scheme II numerically, we use (26), 



(59) 



0, 



7]7T 

T 



arcsm 



Pi 



7X/2' 



77 = 1,1/2,1/4, 
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instead of p8| ), and ([27]) with (a) 77 = 1 as well as with (b) r\ = fj = 1, 1/2, 1/4. Condition ([29]) is 
satisfied if and only if 77 = 77. 

The results for case (a) at 7 = 1 , 20 are shown in Figure [4] and the results for case (b) with 
77 = fj at 7 = 1, 20 are displayed in Figure [5} The slight degradation in performance for e.g., 77 = 
1,77 = 1/4,7 = 1 results from the violation of ( |29| ) (and hence (|23[)) which affects the performance 
significantly for small 7. For large 7, it has been shown in [12J that the sampling angles can be 



chosen independently of the incident angles to maintain a high level of performance (Figure [7]). 
On the other hand, when (29) holds, the performance is essentially independent of both 7 and 77, 
Figure [5} 



For Scheme III we use equally spaced incident angles Q\ G [ir/6,ir/3], (20) and 



(60) 



9i + ei = 2<fn + w, fj = 1,1/2,1/4. 



As shown in Figure [6j the performance is essentially independent of 77. 

We demonstrate numerically the high-frequency SIMO schemes analyzed in [12\ and compare 
their performance with that of the multi-shot SISO schemes presented above. A case in point 



would be to use (26) and (28) but for a fixed incident angle, say 61 = 0, V7 (In this case (29) is 



almost certainly violated). In |12j it is established that the SIMO schemes with sampling angles, 
independent of the incident angles, achieves a high performance in reconstructing a sparse target 
with a sufficiently high-frequency probe wave (i.e. 7 3> 1). The success probabilities of the SIMO 
schemes for 7 = 1, 20, 200 are calculated and plotted in Figure [7| Consistent with the theory |12| . 
the low frequency case with 7 = 1 has the worst performance. Clearly the performance of the 
SIMO schemes improves with 7 and in the limit 7 > 1 approaches that of the multi-shot SISO 
Scheme II. There is negligible difference between the performances for 7 = 20 and 7 = 200 both of 
which follow closely that of the SISO Scheme II with 7 = 1 (black dashed line in Figure [7]). 

Finally we demonstrate the reconstruction with the Littlewood-Paley basis. Since the scattering 



matrix is block-diagonalized by the proposed sampling scheme fl43b-(44) and (46)- (47), we consider 



targets of single pair of dyadic scales. In this simulation, the reconstruction is carried out for a 
target on the scales (1,32) (p = (0,5)) with 42 nonzero coefficients (s p = 42) among 441 possible 



modes (m p = 10, 21 modes in each coordinate) by using 81 samples 



rir 



81). The target and its 



reconstruction in the domain [—10, 10] x [—200,200] are color-coded and displayed in Figure pi 



6. Conclusion and discussion 

We have proposed, analyzed and numerically tested several multi-shot SISO sampling schemes 
which transform the scattering matrix into the random Fourier matrix in the case of point and 
localized extended scatterers and the block-diagonal form of random Fourier matrices in the case 
of distributed extended targets. 

In the case of point scatterers, these sampling schemes are either multi- frequency band limited 
(I) or single frequency outside band (II). For Scheme I , the sampling direction is the backward 
direction analogous to synthetic aperture radar (SAR) while for Scheme II in the high frequency 
limit the optimal sampling direction is in the forward direction, analogous to the X-ray tomography. 



Both schemes produce nearly the same recovery probability with the resolution given by (30), i.e 

£= — 
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Figure 8. Imaging of an extended target of the scales (1, 32): The left plot is the 
exact structure; the middle plot is the the reconstructed profile; the right plot shows 
the (round-off) error. See the color bars for the meaning of the color code. 



We have extended this approach to the case of localized extended targets by interpolating from 
grid points. We have formulated the approximate scheme as inversion with noisy data. In particular, 
we have derived an explicit error bound for the simple piece-wise constant interpolation. 

In the case of distributed extended targets, the block-diagonal form of the scattering matrix in 
the Littlewood-Paley representation means that different dyadic scales of the target are decoupled 
and can be imaged scale-by-scale separately by our method. Moreover, we can determine the 
coefficients in the Littlewood-Paley expansion (40) for scales up to uj/(2it) by using probes of any 
single frequency co. The disadvantage of the Littlewood-Paley basis is that a localized target has 
slowly decaying coefficients and hence is not compressible in this basis. 

The SIMO schemes in which the scattered field of an incident angle is measured at multiple 
sampling angles have been studied in [I2j HH Ej. Except for the special case of the periodic 
scatterers lying on a transverse plane |13| . it is not known if the sensing matrices of the SIMO 
schemes in general satisfy the RIP or not. In this case, the approach based on the notion of 
incoherence is taken to analyze the SIMO schemes \12\ I14j . This approach is generally more 
flexible and should be applicable to the SISO schemes considered here. 

The main advantage of the SIMO schemes is that in the one-shot setting (one incident field) the 
inverse scattering problem can be solved exactly without the Born approximation by inverting an 
auxiliary nonlinear system of equations [12J. We are working to extend the idea to the imaging 
methods with multi-shot measurements. 

On the other hand, the SISO schemes with the RIP tend to have a better performance which 
can be matched by that of the SIMO schemes only at high frequency as demonstrated in Figure 
[7j In the high frequency regime the sampling angles can be chosen independently of the incident 
angles [12]. 

15 



Acknowledgement. I thank my students Hsiao-Chieh Tseng (Figures [3||Tj) and Wenjing Liao 
(Figure [8]) for preparing the figures. 

References 

[1] M. Born and E. Wolf, Principles of Optics, 7-th edition, Cambridge University Press, 1999. 

[2] A.M. Bruckstein, D.L. Donoho and M. Elad, "From sparse solutions of systems of equations to sparse modeling 
of signals," SI AM Rev. 51 (2009), 34-81. 

[3] E. J. Candes, "The restricted isometry property and its implications for compressed sensing," Compte Rendus 
de I'Academie des Sciences, Paris, Serte I. 346 (2008) 589-592. 

[4] E. J. Candes, J. Romberg and T. Tao, "Robust undertainty principles: Exact signal reconstruction from highly 
incomplete frequency information," IEEE Trans. Inform. Theory 52 (2006), 489-509. 

[5] E.J. Candes, J. Romberg and T. Tao, "Stable signal recovery from incomplete and inaccurate measurements," 
Commun. Pure Appl. Math. 59 (2006), 120723. 

[6] E. J. Candes and T. Tao, " Decoding by linear programming," IEEE Trans. Inform. Theory 51 (2005), 42034215. 

[7] E. J. Candes and T. Tao, " Near-optimal signal recovery from random projections: universal encoding strate- 
gies?," IEEE Trans. Inform. Theory 52 (2006), 54-6-5425. 

[8] S.S. Chen, D.L. Donoho and M.A. Saunders, "Atomic decomposition by basis pursuit," SIAM Rev. 43 (2001), 
129-159. 

[9] D. Colton and R. Kress, Inverse Acoustic and Electromagnetic Scattering Theory. 2nd edition, Springer, 1998. 
[10] W. Dai and O. Milenkovic, "Subspace pursuit for compressive sensing: closing the gap between performance and 

complexity," |arXiv:0803.08TT| 
[11] I. Daubechies, Ten Lectures on Wavelets. SIAM, Philadelphia, 1992. 

[12] A. Fannjiang, "Compressive inverse scattering I. High frequency SIMO measurements, " arXiv: 0906.5405 
[13] A. Fannjiang, "Compressive imaging of subwavelength structures," SIAM J. Imag. Set. 2 (2009), 1277-1291. 
[14] A. Fannjiang, P. Yan and Thomas Strohmer, "Compressed remote sensing of sparse objects," |arXiv : 0904 . 3994 
[15] H. Rauhut, "Stability results for random sampling of sparse trigonometric polynomials," preprint, 2008. 
[16] M. Rudelson and R. Vershynin, "On sparse reconstruction from Fourier and Gaussian measurements," Comm. 

Pure Appl. Math. Ill (2008) 1025-1045. 
[17] V. N. Temlyakov, "Greedy approximations," Acta Num. (2008), pp. 235409. 
E-mail address: fannjiang@math.ucdavis.edu 

Department of Mathematics, University of California, Davis, CA 95616-8633 



16 



